Phonon Density-of-States in MgB 2 
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We report inelastic neutron scattering measurements of the phonon density-of-states in Mg 11 B2, 
which has a superconducting transition at 39.2 K. The acoustic phonons extend in energy to 36 meV, 
and there are highly dispersive optic branches peaking at 54, 78, 89 and 97 meV. A simple Born-von 
Karman model is able to reproduce the mode energies, and provides an estimate of the electron- 
phonon coupling of A ~ 0.9. Furthermore, the estimated boron and magnesium contributions to the 
isotope effect are in qualitative agreement with experiment. The data confirm that a conventional 
phonon mechanism, with moderately strong electron-phonon coupling, can explain the observed 
superconductivity. 

PACS numbers: 74.25.Kc, 63.20.Kr, 74.62.-c, 78.70.Nx 



There has already been rapid progress in our under- 
standing of the properties of MgB2, since the first an- 
nouncement of superconductivity with a critical temper- 
ature of 39K § . Although some have speculated that the 
high value of T c must require a novel mechanism for su- 
perconductivity p], initial experimental reports suggest 
that a conventional phonon mechanism is sufficient. Both 
the observation of the isotope effect p,0 and a simple 
BCS form for the energy gap measured by point-contact 
tunneling p] are consistent with phonon-mediated su- 
perconductivity, with the high transition temperature 
made possible by the high frequency of boron vibrational 
modes. Initial estimates of phonon frequencies, from 
electronic structure calculations [^-^), predict that they 
are consistent with the observed critical temperature, al- 
though An et al J5J and Kong et al |8) predict that only 
a few modes have significant electron-phonon coupling. 

The strength and frequency-dependence of the 
electron-phonon coupling is determined by measurements 
of both the electron-phonon spectral density, a 2 (ui)F '(u>) , 
and the bare phonon density-of-states (PDOS), F(u>). 
a 2 (w)F(ti;)can be derived by inversion of tunneling data, 
while F(w) is derived from inelastic neutron scattering 
measurements. In systems where both sets of data are 
available, a(ui) usuallly varies smoothly with frequency 
so the bare phonon density-of-states provides useful esti- 
mates of the frequency moments that determine T c ||. 

We report measurements of the phonon density-of- 
states on an isotopically enriched polycrystalline sam- 
ple of Mg 11 B2 |[o| using time-of-flight neutron spec- 
troscopy at the Intense Pulsed Neutron Source, Argonne 
National Laboratory. They confirm that boron vibra- 
tional frequencies extend as high as 100 meV. The acous- 
tic modes peak at 36 meV, while there are dispersive op- 
tic branches at 54, 78, 89 and 97 meV, in good agreement 
with predictions based on electronic structure calcula- 
tions [^-^]. A simple Born-von Karman (BvK) model, 
which reproduces the energies of the main spectral fea- 
tures, is used to determine the logarithmic frequency mo- 



ment of the PDOS, which provides a preliminary esti- 
mate of the electron-phonon coupling, A ~ 0.9 (assuming 
fx* = 0.1). We also predict the scale of PDOS shifts 
with isotopic mass, providing estimates of the magne- 
sium and boron contributions to the isotope effect, al- 
though the values depend on the definition of the PDOS 
used. Given the uncertainties in this analysis, the agree- 
ment with experiment is satisfactory, and allows us to 
conclude that the superconductivity is consistent with a 
conventional phonon mechanism with moderately strong 
electron-phonon coupling. 

The MgB2 sample was synthesized from a high purity, 
3mm diameter Mg rod and isotopic 11 B (Eagle Picher, 
98.46 atomic % n B). The Mg rod was cut into pieces 
about 4mm long and mixed with the 200 mesh n B pow- 
der. The reaction was done under moderate pressure 
(50 bars) of UHP Argon at 850°C. At this tempera- 
ture the gas-solid reaction was complete in about one 
hour. The sample (approximately 20 g) was contained 
in a 30 ml tantalum crucible with a closely fitting cover. 
There was surface reaction between the Ta crucible and 
the reactants at the synthesis temperature but the prod- 
ucts of this reaction could be mechanically removed. X- 
ray diffraction showed some small, unidentified impurity 
peaks in the resultant powder, but the impurity concen- 
tration is small and should have no effect on the PDOS 
measurements. 

The neutron scattering experiments were performed on 
the Low Resolution Medium Energy Chopper Spectrom- 
eter (LRMECS), using incident energies of 200 meV, 130 
meV, and 50 meV. LRMECS has continuous detector 
coverage with scattering angles from 2.4° to 117.5° so 
that a broad range of wavevector transfers are measured 
simultaneously, e.g. 2 to 12 A -1 at an energy trans- 
fer of 60 meV. This is extremely important in accurate 
determinations of the PDOS when the scattering is pre- 
dominantly coherent, as it is in Mg 11 B2. The coherent 
(incoherent) scattering cross sections are 3.63 barns (0.08 
barns) and 5.56 barns (0.21 barns) for Mg and n B re- 
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FIG. 1. (a) The generalized phonon density-of-states, 
GDOS = J2i=M a ,B Gi ex P (-2Wi)Gi(u)/Mi, determined from 
the inelastic neutron scattering spectrum of MgE>2 measured 
on LRMECS at 8 K with an incident neutron energy 130 
meV integrated over wavevectors from 9 to 12 A -1 . The 
solid line is the estimated multiphonon contribution, (b) 
The calculated GDOS based on the BvK model described 
in the text convoluted with the resolution function of the 
spectrometer. The dashed lines represent the projections of 
the partial contributions to GDOS; the in-plane boron modes 
(long-dashed), out-of- plane boron modes (dotted), in-plane 
magnesium modes (dash-dotted), and out-of-plane magne- 
sium modes (dash-double-dotted). 

spectively. Unless the scattering is purely incoherent, 
it is necessary to average the measured neutron scatter- 
ing over a large volume of reciprocal space in order for 
the resulting measurements to reflect the true phonon 
density-of-states. 

The neutron data taken with an incident energy of 
200 meV shows that the PDOS extends up to 100 meV. 
Above that energy, the scattering only comprises a broad 
multiphonon tail. Figure [l] shows data taken at 130 meV 
at 8K. The peak in the acoustic PDOS occurs at about 
36 meV, while there are broad optic bands at 54, 78, 89 
and 97 meV. The spectrometer resolution decreases with 
energy transfer from 5.9 meV at hui = 40 meV to 3.4 
meV at hu> = 90 meV. The large width of these features 
is not just a resolution effect, but arises from dispersion 



of the optic phonon branches. We have also performed 
measurements at 30K, 50K, and 100K. The only statisti- 
cally significant difference between these temperatures is 
a small suppression of the peak at 54 meV below T c . In 
a polycrystalline average, this is a very small effect, but 
it could still represent a significant change in an individ- 
ual phonon branch. Interestingly, this is the mode that 
An and Pickett predict to have the strongest electron- 
phonon coupling Q], but such conclusions need to await 
single crystal phonon dispersion measurements. 

In a multicomponent system, the inelastic neutron 
scattering law in polycrystalline samples is given by 

S(Q,u>) = ^^exp(-2Wi)^M[nH + l] 

i=Mg,B 1 

(1) 

where Oi and Mi are the neutron scattering cross sec- 
tion and atomic mass of the i th atom and n(uj) = 
[exp (huj/kBT) — 1] _1 is the Bose population factor. The 
generalized PDOS, G(v) = where Gi(uj) is 

defined as 

= jnJ2 i e ^'< - <$] (2) 

which, in turn, defines the Debye- Waller factor, Wj(Q) 
through 

wm = % f°^W) + i] (3) 

2Mi J lu 

u>(j, q) and ei(j, q) are the frequencies and eigenvectors, 
respectively, of the phonon modes. Because the Mg and B 
partial contributions are weighted by cr/M, the measured 
spectrum does not represent the true PDOS, although 
the stronger cross section and lighter mass of the boron 
atoms ensure that they make the dominant contribution. 

We have analyzed the data in terms of a simple BvK 
model, that is sufficent to describe the principal features. 
MgB2 has a hexagonal structure with the lattice parame- 
ters a = 3.082 Aand c = 3.515 A Jn|. The lattice is com- 
posed of a honeycomb of B atoms separated by hexagonal 
Mg layers, leading to comparisons with graphite interca- 
lation compounds (t]] . The nearest-neighbor distances are 
3.082 and 3.515 A for Mg-Mg pairs in the basal plane and 
along the c-axis, respectively, 2.501 A for Mg-B pairs, and 
1.780 and 3.515 A for B-B pairs in the basal plane and 
along c-axis, respectively. 

Assuming central forces, each atom-atom interaction 
can be represented by two parameters, the longitudid- 
inal and transverse force constants. We have included 
only nearest-neighbor interactions between Mg-Mg and 
B-B atoms in the basal plane and along the c-axis, 
and between Mg-B pairs in neighboring layers. The re- 
sults of the calculation, which was performed using the 
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FIG. 2. The bare phonon density-of-states, F(w), of 
Mg 11 B 2 (solid line) and Mg 10 B 2 (dotted line), calculated from 
the BvK model described in the text (normalized to 1.0). 
The inset shows neutron data taken with an incident energy 
of 50 meV showing the split peaks in the acoustic densi- 
ties-of-state. 

UNISOFT program fl^l , are shown in Fig. 1(b), where 
the in-plane and out-of-plane projections of the partial 
contributions to G(u>) from boron and magnesium modes 
are also plotted. 

It is possible to reproduce the frequencies of the peaks 
in the PDOS with only six force constants. These are the 
longitudinal and transverse in-plane B-B interactions (90 
and 7.5 N/m respectively), the longitudinal out-of-plane 
B-B interaction (20 N/m), the longitudinal in-plane and 
out-of-plane Mg-Mg interactions (both 0.6 N/m) and the 
longitudinal and transverse Mg-B interactions (70 and 2 
N/m respectively). The remaining transverse nearest- 
neighbor interactions had a negligible effect on the mode 
frequencies. 

The frequencies of the main peaks in the PDOS are 
well-represented by this model, but it is clear that the 
widths of the peaks are not, i.e. we have underestimated 
the dispersion of the optic branches. It would be pos- 
sible to improve the agreement by the introduction of 
additional parameters representing longer-range interac- 
tions, but the results of such an analysis would be mean- 
ingless at the present time. Complete measurements of 
the phonon dispersions on single crystals would allow the 
model to be refined, but the present analysis is sufficient 
to provide good estimates of the frequency moments of 
the PDOS that determine the phonon contribution to T c . 

Fig. 2 shows the resulting bare PDOS, F(u>) defined 

as 
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6[u-u)(J, q)] 



(4) 



This is similar to equation g but does not include the vi- 
brational eigenvectors. In a Bravais lattice, the two quan- 
tities are equivalent, but in multicomponcnt systems with 



very different atomic masses, G(ui) will tend to weight 
the PDOS to higher frequencies because of the smaller 
atomic displacements of the heavier atoms. 

We now discuss the implications of this data for un- 
derstanding the superconducting transition temperature. 
This may be estimated using the Allen-Dynes equation 



= -yy ex P 



1.04(1 + A) 



A-/i*(l + 0.62A) 



(5) 



where A is the electron-phonon coupling, fj,* is the 
Coulomb repulsion, and ui n is the logarithmic phonon 
average defined in equation 2.16 of ref. j^]. LUi n is an av- 
erage over a 2 (w)F(w), so we need to make the assump- 
tion that a(u>) is approximately energy-independent, in 
which case it cancels from the equation. However, we 
have a choice of two functions to use. It is more con- 
ventional to use the bare PDOS, F(u>), given by equation 
(Q). However, this function fails to take into account the 
size of the atomic displacements of the different phonon 
modes; they are included in G(w). The electron-phonon 
matrix elements contain the product of the mode eigen- 
vectors and the gradient of the crystal potential ||, so 
it is likely that G(lu) represents a better approximation 
to a 2 (uj)F(uj). In fact, the choice of F(u) or G(w) only 
affects significantly the calculation of the isotope effect, 
which we discuss below. 

The Allen-Dynes equation can be inverted to estimate 
the value of A for a given value of T c and /i*. From 
integrating over the model calculation of G(w), we esti- 
mate that LUin is 57.9 meV [53.8 meV when determined 
from F(w)]. Fig. 3 shows calculations of T c using three 
typical values of /x* compared to the measured value of 
Mg 11 B 2 . This shows that the measured value of T c is 
consistent with a phonon mechanism with values of A in 
the range 0.6 to 1.2. A = 0.9 at the most frequently cited 
value of \i* =0.1. u>i n nearly coincides with the energy of 
the Fig mode that An and Pickett predict will dominate 
the electron-phonon coupling 0. It is impossible to de- 
termine, on the basis of our measurements, whether the 
entire phonon spectrum is involved in the superconduc- 
tivity or only a subset that occurs close to the average 
energy. 

Finally, we can use the BvK model to estimate the iso- 
tope effect. The dotted line in Fig. 2 shows a calculation 
of the PDOS in Mg 10 B2, illustrating the significant en- 
ergy shifts in all the high frequency boron modes. Sub- 
stituting 26 Mg for 24 Mg produces much smaller shifts. 
In the absence of strong coupling effects, the total iso- 
tope effect is given by the sum of the partial isotope ef- 
fects, (3 = J2i=Mg,B0i> wli ere $ = -S lnT c /£ lnM* = 
— S lnuJin/6 InMi and should always equal 0.5 if fi* = 
H . More accurate predictions suggest that it can fall 
below this value when A becomes small, but not at the 
value of 0.9 that we estimate. If we use G(ui) to estimate 
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FIG. 3. Variation of the superconducting transition tem- 
perature with A using the Allen-Dynes equation and the value 
of the logarithmic phonon average, u>i n = 57.9 meV, de- 
termined from the generalized PDOS. Three values of the 
Coulomb repulsion, fi* — 0.0 (dashed line), 0.1 (solid line), 
and 0.2 (dotted line), are plotted. The dash-dotted line is the 
measured value of T c in Mg 11 B2. 



Pb and 0Mg, we obtain values of 0.42 and 0.08 respec- 
tively, whereas using F(w), we obtain values of 0.28 and 
0.22, respectively. Reported value of (3b are in the range 
0.28(3) to 0.30(2) Hinks et al have measured an 

average value of /3m s = 0.018(8) Q. The experimental 
values are slightly lower than our calculations but the 
difference may not be significant given the limitations of 
our analysis. 

In conclusion, we have measured the phonon density- 
of-states in Mg 11 B2, and used the data to produce a sim- 
ple BvK force-constant model of the vibrational spectra. 
Our analysis allows us to conclude that the measured val- 
ues of the transition temperature and the isotope effect 
are consistent with a conventional phonon mechanism for 
the superconductivity with moderately strong electron- 
phonon coupling. It will be extremely useful to compare 
our data with tunneling measurements of the electron- 
phonon spectral density, which determine whether the 
whole phonon spectrum is involved in the superconduc- 
tivity, or only a few specific modes. 
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